First principles modeling of tunnel magnetoresistance of Fe/MgO/Fe trilayers 
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By carrying out density functional theory analysis within the Keldysh non-equilibrium Green's 
functional formalism, we have calculated the nonlinear and non-equilibrium quantum transport 
properties of Fe/MgO/Fe trilayer structures as a function of external bias voltage. For well relaxed 
atomic structures of the trilayer, the equilibrium tunnel magnetoresistance ratio (TMR) is found to 
be very large and also fairly stable against small variations in the atomic structure. As a function 
of external bias voltage, the TMR reduces monotonically to zero with a voltage scale of about 
IV, in agreement with experimental observations. We present understanding of the nonequilibrium 
transport properties by investigating microscopic details of the scattering states and the Bloch bands 
of the Fe leads. 
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Since the prediction and elegant physics explanation^, 
2] that magnetic tunnel junction (MTJ) of Fe/MgO/Fe 
trilayer structure may have very high tunnel magnetore- 
sistance (TMR), MgO based MTJ has progressed at a 
rapid pace in recent years and produced the highest mea- 
sured TMR at room temperature: several groups0,0| re- 
ported TMR ratio in the range of 180% to 250%. TMR 
effect presents an excellent opportunity for spintronics, it 
is the key to magnetoresistive random-access- memory [5|, 
programmable logic elements 6], and magnetic sensors. 
The high TMR values in MgO based MTJs have gener- 
ated great excitement for practical applications. 

Atomistic calculations 0, Q have so far played a vital 
role in elucidating the reason behind the observed large 
TMR in Fe/MgO/Fe MTJsfl E H 0. There are, how- 
ever, a number of important issues remain to be under- 
stood from atomic first principles. Most existing work 
predicted^, TMR to be greater than 1000%, experi- 
mental data are still lower. More seriously is that exper- 
imental data on MgO based MTJs show a monotonically 
decreasing TMR as a function of applied bias voltage 0,3 
and it reduces to zero when bias is about one volt. To 
the best of our knowledge, there have been two atom- 
istic calculations of bias dependence of TMR for MgO 
barriers 011], both used the Korringa-Kohn-Rostoker nu- 
merical technique. Ref0 predicted a substantial increase 
of TMR versus bias for the asymmetric system analyzed 
there, while RefJ3 found a roughly constant TMR, a de- 
caying TMR, or an entirely negative TMR versus bias de- 
pending on atomic structures of the interface. The origin 
of these differences were not clear. Earlier theory Q on 
AI2O3 based MTJs has attributed small bias dependence 
of magneto-resistance to magnon scattering. Given the 
extreme importance of MgO based MTJ in near future 
spintronics and the accumulated experimental data, fur- 
ther quantitative understanding on quantum transport in 
Fe/MgO/Fe at finite bias is urgently needed. 

Here we present a first principles atomistic analysis 
of nonlinear and non-equilibrium quantum transport in 



Fe/MgO/Fe MTJ. We use a state-of-the-art quantum 
transport techniaue |l0l ITU which is based on real-space, 
Keldysh nonequilibrium Green's function (NEGF) for- 
malism combined with density functional theory (DFT). 
The basic idea of the NEGF-DFT formalismfiq is to 
calculate device Hamiltonian and electronic structure 
by DFT, populate this electronic structure using NEGF 
which properly takes into account nonequilibrium quan- 
tum statistics, and deal with open device boundaries 
directly using real-space numerical techniques. The 
power of NEGF-DFT methods have already been demon- 
strated by direct quantitative comparison to experimen- 
tal data[I3|- Our results show that for fully relaxed 
atomic structure of the Fe/MgO/Fe device, the equi- 
librium TMR ratio reaches several thousand percent — 
consistent with previous theoretical resultsQ, hj. This 
value is also found to be fairly stable against small varia- 
tions in atomic structure. We found that the TMR ratio 
is monotonically quenched by the bias 14 with a scale of 
about one volt. The microscopic details of these trans- 
port features can be understood by the behavior of bias 
dependent scattering states. 

The MTJ is schematically shown in Fig. la where a 
number of MgO(100) layers is sandwiched by two Fe(100) 
leads. The MTJ is periodic in the x — y direction 
while the leads also extend to z = ±00 (transport di- 
rection). For DFT, wc use standard norm-conserving 
pseudopotentials[Fy and an s, p, d double-zeta LCAO 
basis set^Ej. The exchange-correlation potential is 
treated at the LSDA level. In order to accurately de- 
termine transport properties of the MTJs, we found that 
special care must be given to the pseudopotentials and 
basis sets. In our work, these inputs were carefully con- 
structed to accurately reproduce the electronic structures 
of Fe, MgO, and periodic lattice of the Fe/MgO interface 
obtained by full potential LAPW method^]. The latter 
comparison is shown in Fig. lb. 

For the two probe MTJ simulation, we found that 
20 x 20 (k x ,k y ) points suffice to sample the 2D trans- 
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FIG. 1: (color online) (a) Schematic plot of a two probe 
Fe(100)/MgO(100)/Fe(100) device. The system has infinite 
extent in the (x,y) direction with a lattice constant of 2.82A, 
and extends to ±00 in the z-direction. (b) Band structure 
of a periodic • ■ -Fe/MgO/Fe/MgO- ■ ■ lattice obtained using 
optimized LCAO pseudopotentials and basis sets compared 
to that from full potential LAPW method. A good agreement 
is found to be necessary in order to carry out the NEGF-DFT 
analysis for the two probe Fe/MgO/Fe devices. 



verse Brillouin zone (BZ) for converging the density ma- 
trix on the complex contour energy integration in the 
NEGF-DFT self-consistent analysis [13 • Much denser 
k-sampling of 10 6 (k x ,k y ) points was required to con- 
verge the density matrix for the real energy integra- 
tion of NEGFpiJ, and for computing the transmis- 
sion coefficient by summing over the BZ: T a (E,Vb) — 

J2 k k Ta T ' ky (E, Vb). Here E is the electron energy, 
Vb the external bias. The BZ resolved transmission, 

k k 

To-"' y , is obtained by standard Green's functions tech- 
nique: Tv x ' ky = Tr[Im(E r L )G r Im(i: r R )G a ], where all 
quantities in the trace are functions of transverse mo- 
mentum. Here a =f, J, is the spin index; G r ' a are the 
retarded/advanced Green's function matrices in spin and 
orbital space; and S£ R are the retarded self-energies due 
to the existence of the bulk-3D left/right Fe leads. Fi- 
nally, the spin-current (spin polarized charge current) is 
obtained by I a {V b ) = f dET a (E, V b )[f L (E - fi L ) - 
fn(E — jj, ft)} where /J,l,r is the electrochemical potential 
of the left/right leads. The total charge current is given 
by I = Jf + 1 1 . In our calculations, the atomic structure 
was fully relaxed by the LAPW[16( method between three 
Fe layers on each side of the MgO, with the most remote 
layer of Fe atoms fixed at crystalline positions during re- 
laxation. The x-y lattice constant a of the interface was 
fixed to our calculated one for bcc Fe, a = 2.82A. The 
Fe-0 distance was found to be 2.236Afor a completely 
relaxed structure in agreement with previous studies^}. 

Fig. 2a, b plots the current- voltage (FV) characteristics 
(solid line) for 5-layer MgO device in the parallel mag- 
netization configuration (PC) and the anti-parallel con- 
figuration (APC) of the two Fe leads, respectively. The 
lower insets plot the majority spin-current at small bias 



range. For bias less than 0.8V, the total current remains 
extremely small. At about 1.5V, the device "turns on" 
and current increases rapidly afterward. Such a turn- 
on voltage is consistent with experimental data0, Il7| . 
The spin-currents are shown as the dashed and dotted 
lines for the up- and down-channels (majority-, minority- 
channel) . We found that the initial rise of the current at 
~ 0.8V* in PC is dominated by the down-channel where 
1 1 exceeds 7f by over a factor of eight. This can be ex- 
plained by investigating the transmission coefficients (see 
below). Above ~ 1.5V the spin-currents roughly con- 
tribute equally to the total current. The FV curves for 
a 3-layer MgO device are plotted in the inset of Fig. 2a, b 
and show similar features. 

From the I-V curves we infer a TMR ratio using 
the common optimistic definition: Rtmr. = {Iapc — 
Ipc) I 'I apc , where Iapc, PC ar e the total currents in 
APC and PC respectively. At 14 = when all 
currents vanish, we compute Rtmr. use transmission 
coefficients 0|- From Fig. 2c, for 5-layer MgO device 
Rtmr ~ 3700% at zero bias and declines quickly with 
Vb, essentially vanishes on a scale of about IV. For the 
3-layer MgO we found Rtmr, ~ 850% at zero bias and 
declines to zero with V on a similar bias scale. While the 
experimentally measured TMR has increased dramati- 
cally in the past two years[!| Q, they are still signifi- 
cantly lower than theoretically predicted values here and 
elsewhere 0, 0. It is anticipated that surface roughness 
is playing a major role^. To investigate this effect, we 
generated eight device atomic structures of 5-layer MgO: 
for each device we varied the z-coordinates of the sur- 
face Mg and O atoms from their relaxed positions, by 
a random displacement corresponding to about 1% of 
the bond length. Self-consistent NEGF-DFT analysis is 
carried out for them and the result is shown in Fig. 2d. 
Of these eight atomic structures, the minimum TMR is 
about 3000% while the maximum is ~ 4000%, with an 
average of 3580%. Although the sample size is small, the 
TMR ratio appears rather stable against small interface 
atom displacements. A similar analysis is carried out 
for thirteen 5-layer MgO devices where all the Mg and 
O atoms were displaced randomly by roughly 1% of the 
bond length, the result is in the inset of Fig. 2d. Again, 
even though the sample size is small, the results never- 
theless indicate that small random variations of atomic 
positions in the barrier layer are not sufficient to reduce 
the zero bias TMR to the presently measured experi- 
mental values. Other mechanisms such as oxidization of 
the Fe surface, diffusive impurities and/or defects, are 
likely responsible for experimental TMR values. How- 
ever our results provide a theoretical upper limit which 
does suggest that if a device can be manufactured with a 
high quality interface, it may be possible to achieve even 
higher TMR values than presently knownQ. 

We now investigate noncquilibrium features, namely 
the bias dependence of various transport properties of 
the MTJ. Fig. 2c shows a dramatic quench of TMR by the 
external bias voltage with a scale of about one volt, in 
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FIG. 2: (color online) (a,b) I-V curves for 5-layer PC and 
APC, respectively. Solid line (diamonds): total current; 
dashed line (squares): If, dotted line (circles): If. Upper 
inset: I-V curves for a 3-layer device. Lower inset: majority 
current If vs V& for small ranges of Vj,. (c) TMR vs bias 14 
for a 5-layer device. Inset: TMR for a 3-layer device, (d) 
Histogram of TMR for several 5-layer devices with variation 
in the position of the surface atoms. Inset: Histogram for 
varying all the Mg and O atoms in the device, (e) Transmis- 
sion coefficient T a versus energy E for Vb = 0. E — is the 
Fermi energy of leads. Solid: Tf for PC setup; dashed: Tf 
for PC; dotted: Tf = Tf for APC. Inset: the same transmis- 
sion coefficients at energies between — 3eV and leV. (f) The 
magnitude of the output signal modulation Vout as a function 
of bias. Inset: semi- log plot of zero bias total transmission 
coefficients T versus the number of MgO layers for PC setup, 
indicating an exponential decrease of T. 



agreement with the experimental data|j,|2(j|. The origin 
of the TMR quenching is due to a very fast rise in the 
APC current relative to the PC current as a function of 
bias. We now analyze these features. 

First, the voltage dependence of the total current and 
spin-current (Fig.2a,b) can be understood from the be- 
havior of the transmission coefficient T a . Fig.2e plots 



T a = T a (E) versus electron energy E at zero bias for PC 
and APC of the 5-layer MgO device. In PC, the majority 
carrier transmission Tf (solid line) is smooth and several 
orders of magnitude larger than Tf (dashed line) when 
E is near the Fermi energy of the leads (Ef = 0). By 
analyzing the spin-dependent scattering states [H} of the 
MTJ, we were able to determine which bands of the Fe 
leads contribute to the transmission. We found that Tf is 
dominated by the Ai band of the Fe leads, in agreement 
with Reffl Below —leV, Tf becomes extremely small 
due to the disappearance of the Ai band. The Tf , on 
the other hand, is considerably less smooth because the 
transmission near the Fermi level is mostly dominated by 
interface resonance states |21|. In particular, a large peak 
in Tf appears above E — leV: as E is increased, differ- 
ent Fe bands may participate transport and this peak is 
due to such a contribution. This Tf peak explains the 
much larger minority-channel current than the majority- 
channel current in PC at Vb = 0.8V (Fig. 2a), as already 
noted above. 

Second, for APC, we obtain Tf = Tf for all E at zero 
bias due to the geometrical symmetry of the device (dot- 
ted line in Fig.2e). We found that the BZ resolved total 
transmission, T k *< k *(E,Vb) = T^ x,ky +T J fcx,fcy shown in 
Fig.3c for V b = and Fig.3d for V b = 0.05V, is dominated 
by broad and smooth peaks at around \k x \ = \k y \ = 0.12 
(in unit of 7r /a throughout the paper, where a is the Fe 
lattice constant mentioned above), and there is almost 
no transmission at k x = k y — 0. For V b — 0, Fig. 3c also 
shows that the dominating peaks are surrounded by other 
much sharper peaks. For APC, it is the majority channel 
from one Fe layer traversing the MgO barrier and going 
to the minority channel in the other Fe layer. Figs. 3a, b 
plot the majority and minority electronic band structures 
of Fe near Ef for \k x \ = \k y \ = 0.12, respectively. By pro- 
jecting scattering states with | fc^ | = \ky\ = 0-12 onto the 
Fe bands of Fig. 3a, b, the dominating peaks are found 
to be largely due to channel transmission: they are due 
to the band labelled "1" in Fig. 3a at one Fe contact, 
transmitting to the band labelled "2" on the other Fe 
contact. Our calculations show that this band-to-band 
transmission contributes 2.37 x 10~ 4 to majority channel 
Tf*" y . Other band-to-band transmissions are consid- 

k k 

erably smaller. Similarly, T^ x ' v is mainly contributed 
from band-2 to band-1. Therefore, it is the band-to-band 
transmissions which give almost the entire hight of the 
dominating peaks in Fig. 3c (note Fig. 3c, d plot the total 
BZ resolved transmission) |22j . 

Third, we found that bias voltage has dramatic effects 
for APC. The very sharp peaks in Fig. 3c, which are due 
to interface resonances occurring at zero bias, are com- 
pletely removed by a finite bias of 0.05V, as shown in 
Fig. 3d. Moreover, the dominating peaks become consid- 
erably higher than those in Fig. 3c. We checked that even 
0.01V bias can remove these sharp resonances. Again, 
we found that transmission from band-1 of left Fe lead 

k k 

to band-2 of right lead dominatesT^ m ' y , contributing 
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FIG. 3: (color online) (a,b) Fe bands at \k x \ = \k y \ = 0.12 
versus k z for majority and minority electrons, respectively. 
(c,d) Total BZ resolved transmission coefficient at Ef versus 
k x ,k y , for 5-layer MgO. (c) for V b = 0; (d) for V b = 0.05V. 
The dominant peaks are near \k x \ — \k y \ = 0.12. 

3.2 x 10 -4 to the peaks. This value is considerably larger 
than the value at zero bias, indicating that bias enhances 
the coupling of Fe bands across the MgO barrier for APC, 
causing a fast relative increase in the APC current as a 
function of bias voltage, and is responsible for quenching 
the TMR observed in Fig. 2c. 

An important parameter which has been measured^ 
|2C| is the magnitude of the output signal modulation, 
namely V out = Vb{R.APC ~ Rpc)/Rapc with Vb the ap- 
plied bias, Rapc and Rpc are the junction resistances 
for APC and PC. V ou t puts more weight on information of 
TMR for higher bias, and is plotted as a function of bias 
in Fig.2f. We found that V out increases in a roughly linear 



manner and then bends over at around Vb ~ 0.5 — 0.7V 
where V out is about 350mV. These voltage scales are 
rather similar to the experimental data0, [2£j ■ Finally, 
the inset of Fig.2f shows a semi-log plot of the zero bias 
total transmission at Fermi energy versus four thicknesses 
of the MgO barrier for PC, and the data is in perfect con- 
sistency with the physics of tunneling. 

In summary, we have analyzed non-equilibrium quan- 
tum transport properties of Fe/MgO/Fe MTJs from 
atomic first principles without any phenomenological pa- 
rameter. All the obtained voltage scales for transport 
features are consistent with experimental data, these in- 
clude the turning on voltage for currents, voltage scale 
for TMR quenching, maximum value of V ou t as well as 
the turning over voltage of V ou t- The quench of TMR 
by bias is found to be due to a relatively fast increase of 
channel currents in APC. Very large TMR at zero bias is 
obtained which is stable against small changes of atomic 
positions. The zero bias TMR is, however, expected to 
be sensitive to presence or absence of impurities, oxidiza- 
tion layers and experimental processing procedures. It is 
hopeful that even greater TMR can be obtained experi- 
mentally. 
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